Prevalence of marginally unstable periodic orbits in chaotic billiards 
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The dynamics of chaotic bilhards is significantly influenced by coexisting regions of regular mo- 
tion. Ifere we investigate the prevalence of a different fundamental structure, which is formed by 
marginally unstable periodic orbits and stands apart from the regular regions. We show that these 
structures both exist and strongly influence the dynamics of locally perturbed billiards, which in- 
clude a large class of widely studied systems. We demonstrate the impact of these structures in the 
quantum regime using microwave experiments in annular billiards. 

PACS numbers: 05.45.-a,05.45.Mt 



I. INTRODUCTION 

Chaotic billiards are fundamental paradigms in sta- 
tistical physics and nonlinear dynamics. By connecting 
dynamics with geometry, billiards serve as models to ad- 
dress numerous questions ranging from the foundations 
of the ergodic hypothesis [H, Q and the description of 
shell effects Q to the design of microcavity lasers [3] and 
microwave resonators [a| , among other applications Q . 

A salient feature of billiard systems is that simple ge- 
ometries, such as those in Fig. [1] suffice to give rise to 
a rich variety of dynamical behavior observed in typ- 
ical Hamiltonian systems. But as previously observed 
for specific chaotic billiards, simple geometries may also 
lead to the existence of the so-called bouncing-hall or- 
bits: one-parameter families of periodic orbits exhibiting 
perpendicular motion between parallel walls. Theoretical 
and experimental work on the Sinai [Fig.[T]Jc)] and Buni- 
movich stadium [Fig.[ljd)] billiards have shown that such 
orbits have a major influence on transport properties, 
decay of correlations, and spectral properties H, 0, H, Q. 
This is so because, contrary to the other orbits embedded 
in the chaotic component of the phase space, bouncing- 
ball orbits are only marginally unstable (i.e., perturba- 
tions grow only linearly in time). In general, marginally 
unstable periodic orbits (MUFOs) can be regarded as 
a source of regular behavior that masks strong chaotic 
properties. However, MUFOs are not structurally stable 
and may be destroyed by small changes in the parame- 
ters of the system. Therefore, MUFOs are considered to 
be non-generic and it has long been assumed that they 
could exist only for very special systems, like billiards 
with parallel walls. 

Contrary to this expectation, in this paper we show 
that MUFOs are prevalent in a large class of billiard 
systems. The starting point of our analysis is the ob- 
servation that many of the most widely studied chaotic 
billiards consist of local perturbations of an integrable 
billiard. For concrete examples, consider the chaotic bil- 
liards shown in the right part of Fig. [TJ All these bil- 
liards can be obtained by re-defining the dynamics in the 
gray region of the integrable billiards in Figs, ma)-(c). 



e.g., by introducing a scatterer. It can be shown that 
any orbit (i) lying inside the chaotic component and (ii) 
not interacting with the introduced scatterers will be a 
MUFO. Although bouncing-ball orbits evidently satisfy 
these conditions in the billiards of Figs. [Hd)-(h), the ex- 
istence of such orbits is far from clear in general. Here we 
use geometric and analytical arguments to demonstrate 
the widespread occurrence of MUFOs. Specifically, using 
circular-like billiards as model systems — such as those 
in Figs. [Hf)-(g) — we show that infinitely many fami- 
lies of MUFOs exist for almost all parameter choices of 
the system. We discuss the impact of these structures on 
the dynamics of chaotic orbits as well as the experimen- 
tal observation of MUFOs in the quantum spectrum of 
microwave annular billiards. 




etc... 



FIG. 1: Adding /omi perturbations to integrable billiards, as 
those shown in (a)-(c), one obtains frequently studied chaotic 
billiards, such as those shown in (d)-(h). The gray regions of 
the (a) rectangular, (b) circular, and (c) elliptical billiards are 
defined in such a way that chaotic motion is possible {a < R 
in (b) is the radius of the smallest circle that circumscribes 
all scatterers). MUPOs are shown here to exist in billiards 
such as (d) Sinai [Jl, (e) stadium (f) annular [l^, (g) 
mushroom [ij, and (h) elliptical with scatterers [l5| . 
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The local perturbations described above are typical 
for billiard systems and differ fundamentally from the 
global perturbations considered in smooth Hamiltonian 
systems. In the latter, the KAM theory shows that most 
quasi-periodic orbits of the integrable system survive the 
perturbation, while all periodic orbits with marginal sta- 
bility disappear. Quite the opposite happens in the for- 
mer case: a large set of quasi-periodic orbits disappears 
but there are families of periodic orbits with marginal 
stability that survive the perturbation by "avoiding" in- 
teraction with the localized scatterers. These orbits give 
rise to families of MUPOs detached from regular regions, 
which were previously observed in billiards with paral- 
lel walls Ja, and for specific parameters of the mushroom 
billiard flo, [ll| . Here we consider generic control param- 
eters of a wide class of systems where we characterize the 
MUPOs both theoretically and experimentally. 

The paper is organized as follows. In Sec. [Til we per- 
form a detailed analysis of the existence of MUPOs in the 
annular billiard, a representative example of the class of 
billiards we are interested in. In Sec. lIIII we show the exis- 
tence of an infinite number of different families of MUPOs 
in annular and other circular-like billiards. Our experi- 
mental results on microwave cavities appear in Sec. IIVI 
Finally, our conclusions are summarized in Sec. |Vl 



or form a regular region for |sin(6')| > a, called whis- 
pering gallery. A collision with the scatterer happens 
whenever [l^ 

|sin(6') -(5sin(6'-0)| > r. (1) 

In Fig. [^b this condition is satisfied between the dashed 
lines. In the following we calculate the geometrical con- 
ditions for the existence of MUPOs and we demonstrate 
that typically an infinite number of families q, r/ satisfy 
these conditions. 



(a) (l)=7i/2 




II. ANNULAR BILLIARD 

Annular billiards are defined by two eccentric circles, 
as shown in Fig. [2ta) . For a fixed radius i? = 1 of the 
external circle, the control parameters are the radius r 
and displacement S < 1 — r oi the internal circle, which 
serves as a scatterer. The phase space shown in Fig.ISfb) 
is obtained by plotting the position € [0, 2tt] of the 
collision of the particle with the external circumference 
and the sine of the angle 6 £ [—tt/2, 7r/2] with the normal 
direction right after the collision. In this system, periodic 
orbits of period q and rotation number 77 that collide 
only with the external circumference define star polygons 
of type ((7, 77) , where the integers q and r/ are coprime 
and rj < q/2. A star polygon of type (5,2) is shown in 
Fig. [TJb) and star polygons of types (2, 1) and (5, 1) are 
shown in Fig. [21 Each star polygon belongs to a family 
of orbits of the same type, which is parameterized by (f> 
and has a fixed collision angle sin(0sp) = cos{TTr] / q) . 

For this system, the conditions (i) and {ii) for the exis- 
tence of MUPOs mentioned in Sec. U translate into cross- 
ing the circle of radius a = r + d without colliding with 
the scatterer. Under these conditions, the orbits are em- 
bedded into the chaotic sea but are only marginally un- 
stable (both eigenvalues of the Jacobian matrix equal 1). 
The two orbits shown in Fig. [21 satisfy these conditions 
and hence are MUPOs. We use MUPOs {q, 77) to denote 
the entire one-parameter family of orbits corresponding 
to star polygons {q,rj) that satisfy conditions (i) and (ii). 
Note that these orbits are necessarily periodic because 
non-periodic orbits will either collide with the scatterer 
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FIG. 2: (Color online) Annular billiard for parameters r = 
0.35 and S = 0.5: (a) configuration space and (b) phase space. 
MUPOs correspond to the periodic orbits that cross the circle 
of radius a [dotted line in (a)] but that do not hit the scatterer 
[region between the dashed lines in (b) in which relation ^ 
is satisfied]. The symbols ■ and A indicate, respectively, in- 
dividual orbits belonging to MUPOs (2, 1) and (5, 1). 

Consider MUPOs that encircle the scatterer from out- 
side, like the pentagon-MUPO (5, 1) in Fig. [1 Condi- 
tions for the existence of such outer MUPOs are obtained 
by noting that every star polygon (q, rj) draws an inner 
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regular g-sided polygon, like the pentagon in Fig. [Ijb). 
The radii (d, D) of the inscribed and circumscribed cir- 
cles of this inner polygon are given by c? = cos{TTri/q) and 
D = d/ cos{TT/q). It follows that an orbit of type {q,ri) is 
an outer MUPO {q, rf) if and only if 

cos(7r-^) < cos(7rA) < \ ' \ +r{\ ^^), (2) 

q cos(7r/gj cos(7r/(7j 

where cos(7rA) = a = r + 5. A similar expression is 
obtained for mushroom billiards [l^, ■ 

Inner MUPOs, like the diameter-MUPO (2,1) in 
Fig. [21 exist when 



5 > 



T Tj Tj 1 — Tj 

-— — -— + cos(7r-) + sin(7r-) tan(7r ). 

cos(7r(l - rj)/q) q q q 



Mixed inner-outer MUPOs may also exist for ij > 2. Fam- 
ilies of inner, outer, and mixed MUPOs (5, 2) are illus- 
trated in Fig. [31 

(a) Outer MUPOs (b) Mixed inner-outer MUPOs (c) Inner MUPOs 




outer 



mixed 



FIG. 3: (Color Online) Size w of the families of MUPOs (5, 2) 
in the annular billiard with r = 0.12: (a) outer MUPOs 
for 5 = 0.2, (b) mixed inner-outer MUPO for S = 0.5, and 
(c) inner MUPOs for S = 0.8. All three kinds of MUPOs 
may coexist for a fixed 5. The size w of the families of MU- 
POs in Eq. (|3]) is given by the length of the external arcs. 
MUPOs (5, 2) outside of these regions do not exist since they 
collide with the inner scatterer. 



Unlike the regular regions around stable periodic or- 
bits, MUPOs have zero Lebesgue measure in the phase 
space. The relevant measure is therefore the size of the 
families of MUPOs, given by the length w of the set of 
angles (f> (normalized by 27r) for which an orbit with a 
given {q, rj) exists. In Fig. [2l w is proportional to the 
length of the external arcs of circumference in (a) and to 
the horizontal lines in (b) . For orbits inside the whisper- 
ing gallery one would have w ^ 1. For a given family of 
MUPOs {q, 7]) we calculate w as 



W = Wouter + Wi: 



ed < 1, 



(3) 



where w outer = 1 - qP~ /tt and Winner = qP^/t^ with 
cos(/3=^) = [cos(7r77/g) ± r]/S. For the MUPOs (q, r/) in- 
vestigated in Sec. IIVI below. Wmixed = 0. Figure [31 shows 
the geometrical representation of the terms in Eq. ^ . 



III. INFINITE FAMILIES OF MUPOS 

We now determine the number and values of the dif- 
ferent families of MUPOs {q, 77) 's that exist in a given an- 
nular billiard (r, S). For inner and mixed MUPOs, only a 
finite number of {q, 77) 's exist 16], which can be obtained 
by inspection. On the other hand, we show next that an 
infinite number of outer MUPOs (g, rj) typically accumu- 
late close to the whispering gallery. Let r]{q) denote the 
integer rj for which rj/q — X is minimal and non-negative. 
In the limit g ^ 00 (^^y^ ^ ^) ^ 0+, both inequali- 
ties ([2]) are satisfied if 



A < 



2VT 



J.2 q2 



(4) 



Essentially the same expression is obtained for mushroom 
billiards pJl] and the same scaling on q is expected in the 
case of other circular-like billiards p^ . 

Optimal rational approximants of A = arccos(a) / tt for 
a fixed q are obtained by truncating the continued frac- 
tion representation A = — , 1 = [ai,a2, ■■■], leading 

"l+c,2 + ..- 

to the convergent rj' /q'. The irrational numbers A* for 
which there exists one integer amax such that < Ofmax, 
for all i, are called numbers of constant type. Numbers 
of constant type are difficult to approximate by rational 
numbers and there exist constants Ci , C2 such that [I3 



^2 



< 



q 



< 



Co 



(5) 



for all convergents t]' /q' . Comparing the inequalities (j^ 
and ([5]) we note the same q~^ dependence. Since the 
convergents are the best approximants, the lower bound 
in ([5]) is valid for all rational numbers. Therefore, pro- 
vided that A is a number of constant type, there are 
regions of the control parameters [a7r/(2\/l — r'^) < Ci 
for annular billiards] for which there exist only a finite 
number of families of MUPOs. The numbers of constant 
type are uncountable and dense in the set of real num- 
bers. They have zero Lebesgue measure, however, mean- 
ing that with full probability A belongs to the comple- 
mentary set of irrational numbers for which C2 — > in 
Eq. (O. Therefore, an infinite number of MUPOs exist 
for almost all A and hence for almost parameters (r, 5) . 

The demonstration above can be used in circular- 
like billiards with arbitrary inner scatterers [isl ] to ver- 
ify whether the convergents "q'/q' of A = arccos(a)/7r 
are MUPOs {q',ri') [e.g., satisfy condition ^ in the 
case of annular billiards or Eq. (6) of Ref. [? ] in 
the case of mushrooms]. Typically, an infinite num- 
ber of different families (q, 77) can be found among 
the convergents. For the annular billiard illustrated 
in Fig. [21 for instance, all odd convergents tested are 
MUPOs: (5, 1), (11, 2), (436, 77), (1342, 237), while the 
MUPO (4, 1) is not a convergent. 
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FIG. 4: Experimental length spectrum of the annular billiard 
for 5 — 0.48. Peaks associated to four periodic orbits are indi- 
cated: the shortest one is unstable, the diameter and triangle 
are MUPOs, and the square is inside the regular region. 



IV. EXPERIMENTAL RESULTS 

Having shown that MUPOs are abundant, we now 
study the impact of these structures in quantum exper- 
iments. We use the equivalence between Schrodinger's 
and Helmholtz's equations for flat microwave cavities 
0] to investigate the effect of MUPOs in quantum 
annular billiards. We show that MUPOs are detectable 
and play a prominent role among the periodic orbits. 

A microwave cavity with radius 12.5 cm, height of 

5 mm, and 4 coupling antennas was used in the exper- 
iments. The inner scatterer had a radius of 1.5 cm, 
leading to r = 0.12. The resonance spectra have been 
obtained using a vectorial network analyzer, measur- 
ing the complex amplitude ratio of the input and out- 
put microwave signal of the cavity. For each value of 

6 = 0, 0.08, • • • ,0.88 we measured 10 spectra up to 
10 GHz with a resolution of 100 kHz. Different anten- 
nas and antenna combinations were used to find as many 
resonances as possible. Close lying levels (e.g., split dou- 
blets) were detected as one resonance only due to their fi- 
nite width. However, since the position of those doublets 
can be approximately calculated, a second {not detected) 
eigenvalue could be attributed to the corresponding fre- 
quencies. We justify this procedure by using the high 
precision data obtained with superconducting cavities in 
the experiments described in Ref. [ll|: there the dou- 
blets could be resolved, and we found that the length 
spectrum is stable under small random shifts of one dou- 
blet partner, which allows us to assume doublets to be 
degenerate. Finally, by comparing the number of de- 
tected levels iV(/) below the frequency / to the expected 
number iVvveyi given by Weyl's formula 0], we checked 
that almost all eigenvalues in the considered part of the 
spectrum have been found {N « 150 for each 6). 

Performing a Fourier transform {TT) of the level den- 
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FIG. 5: (Color online) Semiclassical strengths Ssc (lines) com- 
pared to the experimental strengths y (symbols). Ssc is given 
by Eq. ([7]) and expected from periodic orbit theory, while y is 
extracted from the length spectra, as shown in Fig. 2] Three 
orbits are considered: diameter (dotted line and circles), tri- 
angle (solid line and triangles), and square (dot-dashed line 
and squares). The horizontal lines (top) indicate the values 
of 5 for which the corresponding MUPOs exist. 

sity p(/) — '^^Jj:'^'^ we have computed the length spectrum 

\Pfluc{x)\ = |J^r{p(fc) - pWoyl(fc)}| , (6) 

where k = 27r//c. The classical periodic orbits man- 
ifest themselves as peaks located at the corresponding 
orbit length. The length of periodic orbit (g, rj) is given 
by X(^q jjj = 2Rqsm{m^/q). Particularly, for all ^'s we 
consider the peak heights y (strengths) of the diame- 
ter (x(2.i) = 0.5m), triangular (2^(3,1) = 0.63m), and 
square (2^(4,1) = 0.71 m) orbits. The length spectrum for 
5 = 0.48 is shown in Fig. |4l where we indicate addition- 
ally the peak &t x — 0.34 m related to an unstable peri- 
odic orbit. Notice that this peak is much smaller than 
the peaks associated with the MUPOs. 

Using periodic orbit theory, the strength of an orbit in 
a quantum mechanical length spectrum is given by the 
amplitudes of the oscillatory terms in a semiclassical pe- 
riodic orbit summation. We use the trace formula for 
integrable systems to obtain the orbit dependent ampli- 
tudes A = V — — ^''^'^^ , where v — 1 for the diameter 

orbit, and v ^ 2 for all other orbits [3]. The expected 
strength in the case of the MUPOs is 

Ssc{q,ri) ^wA, (7) 

where w is the measure of the entire family, given in 
Eq. ([3]) and illustrated in Fig. [3l In Fig. [5] we compare 
Ssc (lines) with the experimental strengths y (symbols, 
rescaled by a common factor) for different values of 6. 
The dependence of Ssc on 5 is due to the factor w. Over- 
all, the orbits strengths y approximately follow the semi- 
classical behavior Ssc for 5 > 0.3. The deviations can be 
understood qualitatively as the experiment diverges from 
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the semiclassical limit: (1) the finite wavelengths imply a 
spatial uncertainty of the order of the typical width of the 
peaks in the length spectrum; (2) the Fourier transform 
of a finite spectral range generates fluctuations in the 
length spectra (of the order of 10 % of the diameter peak 
height, as seen for x < 0.2 m in Fig. |4|). Nevertheless 
we find that quantum behavior resembles the classical 
behavior in the sense that the data support the use of 
the weighting factors w in the semiclassical strengths in 
Eq. Q. 

V. CONCLUSIONS 

We have demonstrated that MUPOs are prevalent and 
that they must be accounted for in billiard experiments, 
which is a new paradigm that advances previous conclu- 
sions drawn for specific systems [1, H, [13]. In particular, 
MUPOs have not been previously observed in annular 
billiards, despite many theoretical [H, Hi], and exper- 
imental ^1§\ studies, including detailed catalogs of peri- 
odic orbits 23] . We have shown that annular and general 
circular-like billiards typically have an infinite number of 
different families of MUPOs in the chaotic component 
close to the border of the whispering gallery. This should 
be contrasted with the case of billiards with parallel walls 
such as Stadium and Sinai billiards, where only a finite 
number of families of MUPOs exists. 

The above mentioned results can be immediately ex- 
tended to other chaotic billiards defined by local per- 
turbations of integrable systems and are expected to find 
applications in both classical and quantum studies. Clas- 



sically, general arguments on marginal instability can 
be used to show that the resulting stickiness of chaotic 
trajectories to MUPOs generates a universal power law 
pit) ~ t-'^ for the survival probability of nearby parti- 
cles [ll|. This scaling is expected to hold for long times, 
while fluctuations (nonperiodic echoes) occur for short 
times [Tol . lYoj . Studies in the quantum regime have shown 
that orbits with marginal stability are robust to small 
perturbations [1^ and give rise to different transport phe- 
nomena 0]. Recent theoretical studies and microwave 
experiments on chaos assisted tunneling in the annular 
billiard have demonstrated a pronounced effect on the 
tunneling of the so called "beach region" between the 
whispering gallery and the chaotic region [H, |2l[ . Dif- 
ferent mechanisms of dynamical tunneling are currently 
under investigation [24] . and special attention is being 
devoted to mushroom billiards [2^, [2^ . Our results have 
fully characterized the dynamics in the "beach region" of 
annular and mushroom billiards in terms of marginal un- 
stable orbits. The formalization of their contribution to 
dynamical tunneling and a comparison with the existing 
numerical and experimental results are interesting open 
questions. 

Acknowledgments 



We thank M. Brack, B. Dietz, and T. H. Sehgman for 
stimulating discussions and B. Lindner for the careful 
reading of the manuscript. This work was supported by 
DFG (SFB 634), Studienstiftung des Deutschen Volkes, 
and CAPES (Brazil). 



[1] L. A. Bunimovich and Y. Sinai, Commun. Math. Phys. 

78, 479 (1980). 
[2] L. A. Bunimovich, Commun. Math. Phys. 65, 295 (1979). 
[3] M. Brack and R. K. Bhaduri, Semiclassical Physics 

(Addison-Wesley, Reading, Massachusetts, 1997). 
[4] W. Fang, A. Yamilov, and H. Cao, Phys. Rev. A 72, 

023815 (2005). 

[5] A. Richter, in Emerging Applications of Number The- 
ory, The IMA Volumes in Mathematics and its Applica- 
tions Vol. 109, edited by D. A. Hejhal, J. Friedman, M. 
C. Gutzwiller, and A. M. Odlyzko (Springer, New York, 
1999), p. 479. 

[6] H.-J. Stockmann, Quantum Chaos: An Introduction 

(Cambridge University Press, Cambridge, 1999). 
[7] F. Vivaldi, G. Casati, and I. Guarneri, Phys. Rev. Lett. 

51, 727 (1983); P. Gaspard and J. R. Dorfman, Phys. 

Rev. E 52, 3525 (1995). 
[8] A. Backer, R. Schubert, and P. Stifter, J. Phys. A: Math. 

Gen. 30, 6783 (1997). 
[9] H. Primack and U. Smilansky, J. Phys. A: Math. Gen. 

27, 4439 (1994); E. E. Narimanov and A. D. Stone, Phys. 

Rev. Lett. 80 49 (1998). 
[10] E. G. Altmann, A. E. Motter, and H. Kantz, Chaos 15, 

033105 (2005). 



[11] E. G. Altmann, A. E. Motter, and H. Kantz, Phys. Rev. 
E 73, 026207 (2006). 

[12] E. G. Altmann, Intermittent Chaos in Hamiltonian Sys- 
tems, Ph.D. thesis, Wuppertal University, Germany 
(2007). 

[13] N. Saito, H. Hirooka, J. Ford, F. Vivaldi, and G.H. 
Walker, Physica D 5, 273 (1982). 

[14] L. A. Bunimovich, Chaos 11, 802 (2001). 

[15] Our proof of the existence of infinitely many MUPOs 
holds when the scatterers (1) induce chaotic dynamics 
and (2) have a finite number of points at the critical 
caustic [circle of radius a in Fig. [1] (b) in the case of 
circular-like billiards] . 

[16] The width in where such star polygons may exist is 
fixed while the distance between two adjacent vertices is 
2-n/q, i.e., goes to zero with increasing q [see Fig. [21(b)]. 

[17] A. Ya. Khinchin, Continued Fractions (Dover, New York, 
1992). 

[18] C. Dembowski, H.-D. Graf, A. Heine, R. Hofferbert, H. 

Rehfeld, and A. Richter, Phys. Rev. Lett. 84, 867 (2000); 

R. Hofferbert, H. Alt, C. Dembowski, H.-D. Graf, H. L. 

Harney, A. Heine, H. Rehfeld, and A. Richter, Phys. Rev. 

E 71, 046201 (2005). 
[19] B. Dietz, T. Friedrich, M. Miski-Oglu, A. Richter, T. H. 



6 



Seligman, and K. Zapfe Phys. Rev. E 74, 056207 (2006). 

ibtd. 75, 035203(R) (2007). 
[20] F. von Oppen and E. K. Riedel, Phys. Rev. B 48, 9170 

(1993); D. S. Saraga and T. S. Monteiro, Phys Rev. Lett. 

81, 5796 (1998). 
[21] O. Bohigas, D. Bose, R. Egydio de Carvalho, and V. Mar- 

vulle. Nuclear Physics A 560, 197 (1993). S. D. Frischat 

and E. Doron, Phys. Rev. E 57, 1421 (1998); R. W. Robi- 

nett. Am. J. Phys. 67, 67 (1999) 



[22] M. Hentschel and K. Richter, Phys. Rev. E 66, 056207 
(2002). 

[23] G. Gouesbet, S. Meunier-Guttin-Cluzel, and G. Grehan, 

Phys. Rev. E 65, 016212 (2001). 
[24] A. Mouchet, C. Ehschka, and P. Schlagheck, Phys. Rev. 
E 74, 026211 (2006) 



[25] A. H. Barnett and T. Betcke, pre-print nlin/0611059 
[26] A. Backer, in preparation. 



